clear; clc; close all;

% 参数设置
a = 1;
b = 2*pi;
c = 11.5;

% 初始条件
x0 = 0.4;
y0 = 0.8;
z0 = 0.6;

% 迭代次数
N = 100000;

% 预分配数组
x = zeros(1, N);
y = zeros(1, N);
z = zeros(1, N);

x(1) = x0;
y(1) = y0;
z(1) = z0;

% 迭代计算，调用SIMM函数
for i = 1:N-1
    [x(i+1), y(i+1), z(i+1)] = SIMM(x(i), y(i), z(i), a, b, c);
end

% 绘制相图 x-y 平面
figure;
plot(x, y, '.', 'MarkerSize', 1);
xlabel('x');
ylabel('y');
grid on;
